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Abstract 

An acyclic USO on a hypercube is formed by directing its edges in 
sucii as way that the digraph is acychc and each face of the hypercube 
has a unique sink and a unique source. A path to the global sink of an 
acyclic USO can be modeled as pivoting in a unit hypercube of the same 
dimension with an abstract objective function, and vice versa. In such a 
way, Zadeh's 'least entered rule' and other history based pivot rules can 
be applied to the problem of finding the global sink of an acyclic USO. 
In this paper we present some theoretical and empirical results on the 
existence of acyclic USOs for which the various history based pivot rules 
can be made to follow a Hamiltonian path. In particular, we develop 
an algorithm that can enumerate all such paths up to dimension 6 using 
efficient pruning techniques. We show that Zadeh's original rule admits 
Hamiltonian paths up to dimension 9 at least, and prove that most of the 
other rules do not for all dimensions greater than 5. 

1 Introduction 

It is now over 30 years since Khachian showed that linear programming problems 
can be solved in polynomial time |12) . His ellipsoid algorithm and subsequent 
interior point methods are not, however, strongly polynomial time algorithms 
and no such algorithms are known. Pivoting algorithms, such as Dantzig's 
simplex method ^ still offer the possibility of being strongly polynomial. One 
reason for this is that pivoting algorithms follow a path on the graph defined 
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by the skeleton of a polyhedron, and it is widely believed that the diameter of 
this graph is polynomially bounded in the dimensions of the linear program. 

In fact, Hirsch conjectured that the diameter of any rf-dimcnsional polytope 
with n facets, where rt > d > 2, is less than or equal to n — d. Very recently San- 
tos found that this conjecture is false, by exhibiting a polytope with dimension 
d = 43 and n = 86 facets with diameter equal to n — d + 1 [T7]. Nevertheless, the 
belief that the diameter is polynomial is still strong. The subexponential bounds 
of Kalai [TT] and Matousek, Sharir, Welzl [TS] also give further grounds for hope. 
These papers use randomized pivot selection rules, and no deterministic rules 
that achieve these subexponential bounds are known. 

The simplex method is a family of algorithms, with each member of the 
family being determined by a pivot selection rule. In practice, Dantzig's original 
rule works extremely well. However, Klee and Minty [T3] constructed a case 
where the simplex method using this rule follows an exponential length path 
on a family of suitably stretched hypercubes, since called Klee-Minty cubes. 
In fact it visits every vertex, that is, Hamiltonian path, of the hypercube (see 
Figure [T]) . Subsequent research demonstrated that many other pivot rules take 




(a) Dimension 1 or 2 (b) Dimension 3 




(c) Dimension 4 

Figure 1: Klee-Minty cube. 

exponential time on variants of the Klee-Minty examples. Such pivot rules 
include the maximum improvement rule (Jeroslow 10!) and Bland's rule (Avis 
and Chvatal 

There are still some pivot rules for which the behaviour of the simplex 
method is unknown. A particularly interesting set of these rules are the history 
based pivot rules, of which the best known is the least entered rule proposed 
by Zadeh in a 1980 Stanford University Technical Report, that was recently 
reprinted [50]. Very recently, 30 years after it originally appeared, Friedmann 
showed that this rule requires at least sub-exponential time in the worst case 
such as 2*^(^) [7]. A non-trivial upper bound on Zadeh's rule is still unknown. 
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To motivate the least entered rule, Zadeh pointed out a characteristic of 
Klee-Minty examples: some variables pivot very few times and other variables 
pivot an exponential number of times. Zadeh's pivot rule avoids this behaviour 
by making each variable pivot, roughly, the same number of times. In this 
way, it behaves similarly to the random pivot selection rules mentioned above. 
In Section |3] we show that if Zadeh's rule follows a Hamiltonian path on a 
hypercube, then indeed, each variable must pivot an exponential number of 
times in the dimension d of the cube. 

Zadeh's rule differs from the former pivot rules in that it uses information 
from the entire pivot history up to that point. Such pivot rules are called history 
based pivot rules. Besides Zadeh's rule, these include the least-recently basic rule 
d] , the least-recently considered rule [1] , the least-recently entered rule [5] , and 
the least iterations in the basis rule [3] . We remark that for each of these pivot 
rules, there exists no known exponential lower bound. 

In this paper, we study the behaviour of history based pivot rules on an 
abstraction of linear programming known as acyclic unique sink orientations 
(AUSOs) of hypercubes, that were introduced by Szabo and Welzl [TH]. These 
are orientations of the hypercube so that the resulting directed graph is acyclic, 
and each face of each dimension has a unique sink and a unique source, see 
Figure [21 The goal is to find the unique sink of the hypercube. 




Figure 2: 2(a) is a non-USO, 2(b) and 2(d) are AUSO cubes, 2(c) is a USO 



cube which has a cycle. 

As noted in |19j . various optimization problems can be solved using this 
model. The direction on an edge of the hypercube corresponds to an increase 
in the value of an abstract objective function defined on the vertices of the 
hypercube. The concept of abstract objective functions was first introduced by 
Adler and Saigal [1] . The related concept of completely unimodal numberings 
was introduced by Williamson Hoke 0. Although AUSOs need not correspond 
to actual polytopes and objective functions, the notions of linear programming, 
such as bases and pivots, are readily available. Therefore we obtain an abstract 
model on which to observe the behaviour of various history based pivot rules. 

AUSOs on d-dimensional hypercubes have a structure that makes for conve- 
nient notation and terminology. Each vertex is labelled 0, (2^^ — 1) such that 
the binary representation of adjacent vertices' labels differ by exactly one bit. 
Each edge has a direction and an orientation. The direction is given by a num- 
ber 1, d indicating which bit is different between the two endpoints (counted 
right-to- left). The orientation is given by a positive sign (-I-) if the differing bit 
is at the edge's tail and I at its head, and it is given by a negative sign (-) 
otherwise. We will use the terms direction to denote which bit is to be changed 
and signed direction which also specifies the orientation. For emphasis, to spec- 
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ify the direction without sign we use the term unsigned direction (see Figure 
[3]). For some pivot rules only the unsigned direction is important whereas for 




Figure 3: Signed and unsigned direction on a 3-dimensional cube 

others both the orientation and direction are important. 

Although AUSOs do not necessarily correspond to LP digraphs, the above 
vertex labelling can be used to model moving along a path on an AUSO as 
pivoting in the dictionary Xd+i — 1 — Xi for i — 1, ...,d. A pivot +i corresponds 
to a pivot where Xi enters the basis and Xd+i leaves, and a pivot —i corresponds 
to a pivot where Xd+i enters the basis and Xi leaves. This allows the AUSO to 
inherit various pivoting strategies that are defined in terms of LPs. 

The Klee-Minty examples mentioned above can be modeled as AUSOs. In 
fact, they show that Dantzig's original pivot rule for the simplex method leads to 
a Hamiltonian path on an associated AUSO for each dimension d. As mentioned, 
similar results have been found for other pivot rules. In this paper, we investigate 
whether history based pivot rules can lead to Hamiltonian paths on AUSOs. 

Our focus on Hamiltonian paths has the following motivations. Firstly, if 
such paths exist for a given pivot rule they are obviously the worst case ex- 
amples. Secondly, the number of AUSOs is extremely large. Stickney showed 
their are 19 in 3 dimensions [18], Moriyama's program showed there are 12640 
in 4 dimensions [TB], and Matousek [T3] has shown that there are at least 2^ 
AUSOs in d-dimensions. So just listing their degree sequences when d ~ 6 re- 
quires at least 2"^^ steps. Except for extremely low dimensions, it is therefore not 
possible to construct all acyclic USOs. However searching for all acyclic USOs 
which contain a Hamiltonian path greatly reduces the search space. This is due 
to a remarkable indegree characterization due to Williamson-Hoke discussed in 
Section [3] We are able to exploit this property 'on the fly' to eliminate early 
prefixes of Hamiltonian paths that cannot be completed to an acyclic USO. This 
is because the final indegree of each vertex is known as soon as it enters the 
path. The enumeration enabled us to see that in fact most rules do not follow 
Hamiltonian paths, a fact we were then able to prove. Of course proving that 
a pivot rule cannot follow a Hamiltonian path does not say anything about the 
existence or not of other exponential length paths. However searching for these 
is likely to be significantly more difficult. 

The paper is structured as follows. In the next section we define various 
history based pivot rules that have appeared in the literature: Zadeh's original 
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rule, least-used direction rule, least-recently considered rule, least-recently basic 
rule, least-recently entered rule, and least iterations in the basis rule. We also 
give an example that shows they are all different. In Section [3] we develop an 
algorithm that generates all Hamiltonian paths, if any, followed by these history 
based pivot rules. We also provide computational results that show that most 
of these rules do not in fact produce Hamiltonian paths for dimensions up to 
7, except in very low dimensions. In Section U we prove this fact holds for all 
higher dimensions for four of the history based rules we have presented. 

2 History based pivot rules 

In this section we review a number of history based pivot rules that have ap- 
peared in the literature starting with Zadeh's original rule. We also present an 
example to show that the rules all behave differently. The difference of these 
rules can be seen from the difference of the history array h. This array is indexed 
by the 2d directions (or sometimes, all d unsigned directions), and represents 
current historical information required for the given rule. 

Zadeh noticed that the Klee-Minty construction (see Figure[T]) greatly favours 
some directions over others, and designed a new pivot rule to defeat this. 

Zadeh's rule (a.k.a. the least entered rule) |20) : For the entering variable, 
select the improving variable that has entered the basis least often thus far. 
(Figure E]) The history array h is defined on all 2c? directions, and h{t) is 
the number of times the direction t is used. 

In Zadeh's rule, as in others that we will study, there may be ties in selecting 
the entering variable. We will assume that ties may be broken arbitrarily in this 
paper. Note that Zadeh's rule chooses between all 2d variables {d decision vari- 
ables and d slack variables) whereas the next history-based rule chooses between 
the d pairs of decision and slack variables, {xi,Xd+i), each of which defines a 
direction. Directions are not a very useful concept in arbitrary polytopes, as 
no two edges may be parallel, but they are a natural feature of hypercubes and 
are inherited by zonotopes, which are projections of hypercubes. They directly 
inherit the d directions of the hypercube, some of which may no longer appear. 

Least-used direction rule (LUD) For the entering variable, select the 
improving variable whose unsigned direction has been used least often thus 
far. (Figure [5]) The history array h is defined on all d unsigned directions, 
and h(t) is the number of times the direction t is used. 

We now give some other history-based rules that have appeared in the liter- 
ature. We show the paths generated by these rules on the previous example in 
the Appendix. 

• Least-recently considered rule [T : Fix an ordering of the variables vi,V2, 
V2d and let the previous entering variable be Vi. For the entering vari- 
able, select the improving variable that first appears in the sequence 
Wi+i, Ui+2, W2d, wi, Wi-i (or vi,...,V2d if this is the first pivot). The 
history array h is defined on all 2c? directions and is initialized by set- 
ting h(t) to be the rank of t in the given fixed ordering. If direction s 



5 



Figure 4: Zadeh's rule 
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is chosen the array is updated as h{t) -(r- {h{t) — h{s) — 1 mod 2d) + 1. 
The Appendix shows the example of the case when initial sequence is 
{+2, -4, +1, -3, -2, +3, -1, +4} (FigureH]) 

• Least-recently basic rule [Johnson in [4 ]: For the entering variable, select 
the improving variable that left the basis least-recently. The history array 
h is defined on all 2d directions: h{t) is the step number the \t\-th bit of 
the vertex was last 1 if t is positive or was last if i is negative. (Figure 

HH) 

• Least-recently entered rule (a.k.a. least-recently used) |6j: For the entering 
variable, select the improving variable that entered the basis least-recently 
thus far. The history array h is defined on all 2d directions: h{t) is the 
step number when the |t|-th bit of the vertex last changes from to 1 if t 
is positive or from 1 to if t is negative. (Figure [T5|) 

• Least iterations in the basis rule [3 : For the entering variable, select the 
improving variable that has been in the basis for the least number of 
iterations. (Figure The history array h is defined on all 2d directions, 
and is the number of times the |t|-th bit of the vertex is 1 if t is positive 
or if t is negative. 

Note that all the examples illustrate distinct paths on the same AUSO cube. 
In the following section we will describe an algorithm to determine if there 
are any AUSOs that admit Hamiltonian paths for the history based methods 



6 




Figure 5: Least Used Direction rule 
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described in this section. 

3 Searching for Hamiltonian paths on AUSOs 
that follow history based pivot rules 

We developed an algorithm for determining if the various history based rules 
can be made to follow a Hamiltonian path on an AUSO. As noted in the intro- 
duction, it is known that the number of AUSOs is a doubly exponential, so a 
direct search quickly becomes infeasible. We use the fact that we are looking 
for AUSOs with Hamiltonian paths, which greatly reduces the search space. 

3.1 Preliminaries 

Our basic approach is to generate Hamiltonian paths starting with an unoriented 
hypercube, rather than first orienting the cube and checking if it is Hamiltonian. 
Suppose that a cube has a Hamiltonian path labelled with vertices vi,...,vn- 
Then acyclicity implies immediately that each edge of the hypercube ViVj with 
i < j must be directed from vt to Vj. Therefore, given a Hamiltonian path on the 
cube, we can easily construct the unique acyclic orientations for all edges of the 
cube. It still remains to test whether this orientation is an AUSO. Fortunately 
there is an efficient way to do this based on Williamson Hoke's theorem [B]: 
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Theorem 3.1. If an orientation on a d-dimensional cube is acyclic, the following 
conditions are equivalent. 

• The orientation is a unique sink orientation. 

• For k = 0, d there are exactly (^) vertices with indegree k (and hence 
(^) vertices with outdegree k). 

This makes it very easy to check if a given Hamiltonian path appears in an 
AUSO cube: we need only test the degree sequence. Furthermore, we can even 
use this test as the Hamiltonian path is being constructed. Note that when a 
vertex is added to the path its indegree and out-degree are known. Also partial 
degree information is known for unexplored vertices. Therefore if Williamson 
Hoke's condition is violated, we need not complete the construction of the given 
path. This leads to an efficient pruning technique. We also have the following 
interesting corollary. 

Corollary 3.1. In a Hamiltonian path on a AUSO c?-cube starting from vertex 
0, the indegree of a vertex is 1 if and only it is reached by a positive direction 
(or, unsigned direction) that is being used for the first time. 

Proof. Suppose {v,v') is an edge on the Hamiltonian path that uses the direction 
+t for the first time. Then all previous vertices in the path must have zero on 
the t-th bit. However all neighbours of v' on the hypercube except v have one 
on this bit, so they cannot have been visited yet. Therefore the indegree of v' is 
one. Since there are d directions, this yields d vertices on the path with indegree 
one. By Williamson Hoke's theorem this is the entire set of such vertices. □ 

As mentioned in the introduction, Zadeh's rule encourages each variable to 
be used as a pivot variable roughly the same number of times. We make this 
precise in the following result. 

Theorem 3.2. Assume that there is a Hamiltonian path P that follows Zadeh's 
rule on an AUSO of an n-cube. The least-used signed direction is used at least 
— I times. 

n 2 

Proof. We may assume that P starts at vertex zero. Let —t to be the least-used 
direction, and k be the number of times signed direction —t is used. Partition 
the d-dimensional AUSO into two {d — l)-dimensional hypercubes Ci and C2 
where the direction t separates the two (see Figure ^ and P starts in Ci . 



Let rrii, i = 1, ...,n he the number of times that signed direction +i is used in P 
and mn+i,i = 1, ...,n be the number of times that signed direction —i is used 
in P. Since all 2" vertices are visited, we have 

2n 

^ m, = 2" - 1 

1=1 

We know that the minimum value rrin+t — k and rrit — k + I. We can estimate 
the sum in another way by computing as P is followed. Suppose we are at 
vertex v in Ci and follow a signed direction +i with > k + 2. The signed 
direction +t would have been a preferred choice since rrit < k + 1. If signed 
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Figure 6: The two {d— l)-diniensional cubes Ci and C2 separated by the direc- 
tion t. 



direction -\-t was not taken, then it must be that its neighbour in C2 was already 
visited. We call this a blocked pair. A similar analysis holds if v is in C2 and 
a signed direction —i is chosen with m„+i > fc -I- 1. There can be at most 2"~^ 
blocked pairs. So in computing the sum of the rrii along P we have at most a 
contribution ofn(fc + 2) + 7T,(fc + l) — 1 for the unblocked pivots and a contribution 
of at most 2"^^ for the blocked pivots. Therefore 



2n 

E 

i=l 



m, < n{2k + 3) - 1 + 2"~ 



Combining the two expressions for the sum, the theorem follows. 



□ 



Unfortunately Theorem 13.21 onlv holds when the path is Hamiltonian. It is 
possible for a non- Hamiltonian exponential length path to use a signed direction 
as few as zero times! An example is shown in Figure [T] Here we assume that 
Ci and C2 are copies of an AUSO cube with a long path. The resulting cube C 
is easily seen to be an AUSO. Note that since the path is non-Hamiltonian in 
C, vertices unvisited by the path in Ci may be directed into C2. 





*- 



V ir V 





Figure 7: An example of a non-Hamiltonian exponential path where one signed 
direction is never used 



3.2 The algorithm and its vahdity 

In this subsection we describe an algorithm that can generate, up to equivalence, 
all Hamiltonian paths on AUSOs using any of the history based pivot rules 
described in Section [5J In this paper, when we say two paths are equivalent, it 
means they are equivalent up to permutation. In other words, when two paths 



P and Q are equivalent, there is a permutation of cordinates / : {1, 2, • • • , d} ^ 
{1,2,- •• ,d} such that f{P) = Q. Algorithm [T] gives the pseudocode of our 
algorithm. We assume that the Hamiltonian path starts from the vertex labelled 
0. We denote the indegree of the vertex x by indeg{x). For t = ±1, ±d the 



Algorithm 1 Enumerate HP on AUSO-cube with history based pivot rule 



1: path {0}. 

2: if current path is Hamiltonian path then 

3: if path is USO then 

4: output the result 

5; end if 

6; else 

7: m ^ minjg|_|_-^ ^ jg feasible 't^^'') \'n^ove{path.end,t) is not visited } 

8: for all t such that h{t) = m do 

9: if h{t) = and 3t' < t s.t. h{t') = then 

10: continue 

11: else 

12: V <r~ move(j)ath.end,t) 

13: if h{t) ^ and indeg{v) = then 

14: continue 

15: end if 

16: path <r- path + v 

17: renew h 

18: continue searching (from line 2) 

19: recover h 

20: delete path. end 

21: end if 

22: end for 

23: end if 



function move{x, t) returns the neighbour of x using the signed direction t, that 
is, the vertex x + sign{t)2^^^~^ . Note that we focus on this function only when t 
is a feasible move. The array h denotes the history information of the path and 
depends on the pivoting rule. For example, in the case of Zadeh's rule, h{t) is 
the number of times the signed direction t is taken. We claim that the algorithm 
outputs, up to equivalence, all required Hamiltonian paths and that there are no 
duplications. First of all we show that each of the required Hamiltonian paths 
are equivalent to one of the paths output by the program. Below, by 'history 
based pivot rule' we refer to any of the rules described in Section [21 

Lemma 3.1. For every Hamiltonian path P on a c?-cube which can be followed 
by a history based pivot rule, there is a labelling of the cube such that P 
begins with vertex and the order of positive directions first used in P is 
{l,2,...,d-l,4. 

Proof of Lemma \3.1\ We show how to embed P on a d-cube so that it has the 
required properties. We label the first vertex in P as and the initial edge of P 
as direction +1. Continuing, for i = 2, ...,d we consider the first edge of P that 
leaves a face of the cube of dimension i — We define the direction used by this 
edge as +i. This induces a labelling of the cube with the desired properties. □ 
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We remark that all paths produced by Algorithm 1 satisfy the conditions 
of Lemma [3. II due to lines 8-10. Next we will prove that Algorithm 1 does not 
produce duplicate paths. 

Lemma 3.2. Let P and Q be two Hamiltonian paths produced by Algorithm 
1. If there is a bijection (permutation of the cordinates) / : {1,2,- ■■ ,d} 
{1, 2, • ■ • , d} such that f{P) ~ Q then it is the identity mapping, i.e. P = Q. 

Proof of Lemma \3.2\ As remarked, both P and Q satisfy the conditions of 
Lemma 13.11 Since f{P) = Q, both paths must use the fc-th positive direc- 
tion for the first time at the same time. By the lemma this must be direction 
-t-fc, hence / is the identity mapping. □ 

As a consequence of these two lemmas we have the following result. 

Theorem 3.3. Algorithm 1 provides a complete duplicate free list of Hamilto- 
nian paths on AUSO-cubes that follow a given history based pivot rule. 

3.3 Computational results 

We implemented the algorithm and ran it on an Optcron computer with 2. 2 GHz 
CPU, 4x4= 16 processors and 132GB of memory. We were able to do a 
complete enumeration up to dimension 6 and the results are shown in Table [T] 
The other rules refer to the least-recently considered, least-recently basic and 
least iterations in the basis rules. We see that the number of Hamiltonian path 
increases exponentially with Zadeh's least entered rule, whereas it becomes zero 
with the other pivot rules. On the basis of these results we conjecture that, 
except for the least entered rule, such Hamiltonian paths do not exist in any 
dimension greater than 6. We present proofs of these conjectures in the next 
section for all rules except for the least-used direction rule. 

Table 1: The number of Hamiltonian paths produced by history based pivot 
rules 



Dimension 


2 


3 


4 


5 


6 


7 


Least entered rule 


1 


2 


17 


1,072 


3,262,342 


> 10^" 


Least-used direction 


1 


1 


1 


2 








Least recently entered 


1 


1 


1 











Least-recently considered rule 


1 


3 


13 











Least-recently basic rule 


1 

















Least iterations in basis rule 


1 


















We also conducted an experiment to check whether these paths satisfy the 
Holt-Klee condition [9], a necessary condition for realizability of an LP-cube 
which states that every d-dimensional faces have at least d disjoint paths from 
a unique source to a unique sink (See Table [2]). For dimension 7 with least 
entered rule, we have not found any such paths, but the computation was not 
completed due to the long running time. 
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Table 2: The number of Hamiltonian paths produced by history based pivot 
rules which satisfy Holt-Klee condition 



Dimension 


2 


3 


4 


5 


6 


7 


Least entered rule 


1 


2 


12 


79 


360 


? 


Least-used direction 


1 


1 


1 











Least recently entered 


1 


1 


1 











Least-recently considered rule 


1 


3 


12 











Least-recently basic rule 


1 

















Least iterations in basis rule 


1 


















4 Non-existence of Hamiltonian paths 

In this section we prove that, except for the least entered rule and least-used 
direction rule, there are no Hamiltonian paths for the history based pivot rules 
considered except for those shown in Table [TJ 

Theorem 4.1. The least iterations in basis rule and the least-recently basic 
rule do not have any Hamiltonian paths on a d-cube for d > 3. 

Proof. Suppose there is such a Hamiltonian path P for some d-cube. From 
Lemma 13.11 we can verify that the first d edges of P must take the directions 
1, 2, d. Therefore P begins with d + 1 vertices 0, 1, 3, 2'' — 1. The origin 
has indegree and the other d vertices have indegree 1 in the AUSO induced by 
P. At this point, for each of the pivot rules, the direction —1 has the minimal 
value of h among the all the outgoing directions, so the first d + 1 steps have 
to be 0, 1, ...,2'* — 1,2^^ — 2. The vertex 2'* — 2 also has indegree 1, since it is 
not adjacent to any of the other vertices already on P. There are d -I- 1 vertices 
which have indegree 1 in total: 1,3,7, ■• • ,2'' — 1, and 2*^ — 2. This violates 
Williamson- Hoke's condition given in Theorem 13.11 which allows only d vertices 
to have indegree 1. □ 

We remark that this proof also can be used to show that the least-recently 
considered rule cannot have a Hamiltonian path for any d > 3 if the ordering 
begins with +1, +2, • • • , +d, —1. 

Theorem 4.2. The least-recently entered rule does not have any Hamiltonian 
paths on a d-cube for d > 5. 

Proof. Suppose P is a Hamiltonian path produced by Algorithm 1 for the least- 
recently entered rule when d > 5. We will show that P must begin with the 
sequence of vertices Q = {Qi, (52, Qs, Q4} where Qi = {0, 1, 3, • • • , 2''— 1}, Q2 = 
{2'^-l-2'^-^, 2''-l-2^-2-2'^-3, . . . , 2''-i}, Q3 = {2'^-^+2, 2, 6, 14, • ■ • , 2''-2} 
andQ4 = [2'^ ~2-2'^-^,2'^ -2-2"^-^ -2"^-^ ■ ■ ,2'^-^ + 2 + A + S,2'^-^ + 2 + 4,}. 

Q includes the vertices {2'^"i-|-2, 2''"i+2+4+8, 2''"i+2+4} as a subsequence 
and docs not contain the vertex 2"^^^ + 2 + 8. These four vertices lie on a 2-face 
which has two sources, 2'*"^ + 2 and 2"*^^ + 2 + 4 + 8, a contradiction (See Figure 
[5]). It remains to show that P begins as specified. 
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Figure 8: 2-dimensional face with two sources and two sinks 

• Qi = 0, 1, 3 • • ■ , 2^^ - 1. This fonows from Lemma 1511 

• Q2 = S'^ - 1, 2'^ - 2''-2 _ 1^ - 2'^-2 _ 2'*-3 - 1, • • • , 2'^-!. 

We prove this by mathematical induction. For the basic step, we will 
show only 2** — 2''^^ — 1 can come right after 2'' — 1. When we visited 
the vertex 2*^ — 1, all of the bits are 1. It means the next vertex can be 
represented as 2'' - 2*= - 1 = J^tZo 2* - 2*^ (d - 1 > fc > 0). By Corollary 
13.11 vertex J2'i=o ~ ^'^ should have two visited neighbours, one of which 
is obviously the vertex 2'' — 1. In other words, there exists j ^ k such that 

Eto 2^-2'-2^ e {0,1,3- •• ,2'i-l} = [v\3ls.t.v = Y!,=„'^}^{^}- 
Since d > 3 forces J^tZl 2' - 2*= - V not to be equal to 0, Y!lZa 2* " 2'' - 2^ 
should be represented as E'=o 2' Y^tZa 2* — YltZi+i 2' for certain 
Therefore, the (A:, j) equal [d — l,d — 2) or {d — 2,d ~ 1), and d — 1 > k 
requires k — d — 2. (See Figure |9] for the binary representation). 



000...000 — »ooo...ooi — »ooo...oii — >- — » 01...111 

j-i-thbit I 
1...101...1 •< — 11. ..Ill 

Indegree of this vertex 
is more than one, 
which requires i - d-2 



Figure 9: Binary representation for the basic step of Q2 

We can prove the inductive step similarly. If the path is continued by 
2'* - 1, 2"* - 1 - 2'^-^, • • • , 2'^ - 1 - \j2tZd-2-k 2'}) tlie next vertex should 
be equal to ^to 2^-Etl2-fe 2' + 2-'' {d-2-k < j < d-2) or ^to 2'- 
Ef=<i-2-fe 2' -2Uj = d-loi- j < d-2-k). By CoroUaryETl two neigh- 
bours of it arc in |o, 1, 3 • ■ • , 2'' - 1, 2'' - 1 - 2^^'^, • • • , 2'' - 1 - Y.'lZl-2-k 2'} 
Using binary numbers, 2'^-|x;Lo 2''-(2+'=)|-l can be denoted 100 .. . 0011 ... 11, 
where we have fc + 1 Os. (See Figure fTOl for the binary representation). 

• Q3 ^la-^-i + 2,2,6,14,- •• ,2^^- 2} 
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000...000 — » 000...001 — * 000...011 — *■ — » 01...111 

One of them should be flipped I 
10...011...11 -c — — 1011...11 -c — 11...111 

1 

10...011...101...11 

Indegree of this vertex is more 
than one, which requires there 
is no 1-bits between 0-bits. 



Figure 10: Binary representation for the inductive step of Q2 

At the vertex 2'^~^, the history array becomes 

{d + a; (if a; > ) 
l{iix = -d) 
2d+l-x(ifd-l<a;<0) 
(See Table [3]). Ahhough its minimum value is 1, when x — —d, and 



Table 3: The history information of least-recently entered rule. 



Vertex 
(binary) 


direction 


Comment 


+ 1 


-1 


+2 


-2 


+3 


-3 




+d- 1 


-(d-1) 


+d 


-d 


0(000. ..000) 





1 


(1 


1 





1 







1 





1 


initial state 


1(000. ..001) 


2 


1 


(1 


1 





1 







1 





1 




3(000. ..Oil) 


2 


1 


3 


1 





1 







1 





1 






























2"-' - 1(011. ..Ill) 


2 


1 


3 


1 


1 


1 




d 


1 





1 




2" - 1(111. ..Ill) 


2 


1 


3 


1 


1 


1 




d 


1 


d+ 1 


1 


end of Qi 


2" - 1 - 2''-^(l(>l...lll) 


2 


1 


3 


1 


1 


1 




d 


ii+ 2 


d+ 1 


1 






























2"-' + 1(100. ..001) 


2 


1 


3 


2d - 1 


1 


2d -2 




d 


d + 2 


d+ 1 


1 




2''-i(100...000) 


2 


•Id 


3 


2d - 1 


4 


2d -2 




d 


d + 2 


d+ 1 


1 


end of Q2 


2"-' +2(100. ..010) 


2 


2d 


2ii+ 1 


2d - 1 


4 


2d -2 




d 


d + 2 


d+ 1 


1 




2(000. ..010) 


2 


2d 


2ii+ 1 


2d - 1 


4 


2d -2 




d 


d + 2 


d+ 1 


2d+ 2 




6(000. ..110) 


2 


2d 


2d+l 


2d - 1 


2d +3 


2d -2 




d 


d + 2 


d+1 


2d+ 2 






























2"-' - 2(011. ..110) 


2 


2d 


2d+l 


2d - 1 


2d + 3 


2d -2 




3d-l 


d + 2 


d+1 


2d + 2 




2" - 2(111. ..110) 


2 


2d 


2d+ 1 


2d - 1 


2d+ 3 


2d -2 




3d-l 


d + 2 


3d 


2d + 2 


end of Q;i 



the second smallest value is 2, when x = +1, we can not use either the 
direction — d or +1, since they lead to visited vertices. That leads us to use 
the direction +2, whose value is third smallest. The vertex 2''^^ +2 enables 
us to use the direction —d at last. Afterward, to avoid visiting an already 
visited vertex, we have to follow the sequence {2'^~^ + 2, 2, 6, 14, • • • , 2'' — 2} 

• Q4 {2''-2-2'^-^2''-2-2''-2_2<i-3^... ^2'^-ij^2 + A + S,2'^-^ + 2 + A} 

This follows the same reasoning as Qs, that is, using the smallest direction 
which reaches unvisited vertex fixes Qa. Note that direction -|-1 can not 
be used because the destination has already been visited in Q2. 

□ 

For the least entered rule, we can prove the following feature concerning the 
beginning of any Hamiltonian path. 
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Theorem 4.3. For every d dimensional Hamiltonian path using the least en- 
tered rule, different signed directions are used for the first 2d — 1 steps. 

Proof. Let w be a vertex visited during the first 2d — I steps. It is enough to 
show there is a direction t and an unvisited vertex move{v, t) for which h{t) = 0. 

If t > and h[t) — then vertex move{v, t) cannot have been visited yet, so t 
is a candidate direction. Otherwise, if for each i > we have h{t) — 1 and there 
exists at least two negative directions — ti, —t2 such that h{—ti) — ^2) = 0. 
Assume that the direction +ti was used earlier than +^2- So for all vertices 
visited so far, it is impossible to have both the t2-th bit at 1 and the ii-th bit 
at 0. This means that the vertex 'move{v, —ti) has not been visited and —ti is 
a candidate direction. □ 

5 Discussion 

From our computational experiments, Zadeh's least entered rule seems very 
likely to have Hamiltonian paths on AUSO cubes. Using our program, we could 
verify such paths exist up to dimension 9, but did not yet find any for dimension 
10. Furthermore, we could not find any general construction, so this is an open 
problem. Even if such Hamiltonian paths exist, it is not clear whether or not 
they could be obtained on AUSOs that are realizable as polytopes. 

Although we showed that a number of history based pivot rules do not 
admit Hamiltonian paths in general, they may still admit exponential length 
paths. Since our program makes heavy use of the fact that we are searching for 
Hamiltonian paths, we were not able to use it to check this for low dimensions. 

Our computer result for Zadeh's rule allow ties to be broken arbitrarily, as 
does the theoretical lower bound obtained in [7]. It would be interesting to see 
the effects of various deterministic tie breaking rules on these results. 
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7 Appendix 




f9) ^ ^1 jjl 



13) ^ ^ 





Vertex 
(binary) 


direction 


Outgoing directions 
(bold for cliosen) 


+1 


-1 


+2 


-2 


+3 


-3 


+4 


-4 


0(0000) 


3 


1 


1 


5 


6 


4 


8 


2 


+U2,+3:+4 


2(0010) 


2 


6 


8 


4 


5 


3 


7 


1 




6(0110) 


5 


1 


3 


1 


8 


6 


2 


4 




14(1110) 


3 


1 


1 


5 


6 


4 


8 


2 


-3 


10(1010) 


7 


3 


5 


1 


2 


8 


4 


6 


+L-2 


8(1000) 


6 


2 


4 


8 


1 


7 


3 


5 





Figure 11: Least recently considered 
rule 




Figure 12: Least recently basic rule 



Step 
Dumber 


Vertex 

(binai)') 


itectioii 


Ootgoioglrectioiis 
(bold for choseo) 


11 


-1 


+2 


-2 


+3 


-3 


+4 


-4 


1 


i(iiii) 








1 


i 


1 




1 


|I+2,+3,|4 


2 


1(11) 


2 




i 


2 


i 


2 


i 


2 


+2,|U4 


•] 


5(il) 


3 






3 


3 


2 


i 


3 


-1;+2,|4 


4 


13(1101) 


4 




i 


4 


4 


2 


4 


3 


+2,-3 


5 


15(1111) 


5 




5 


4 


5 


2 


5 


3 


-1.-3,-4 


6 


14(1110) 


5 


6 


6 


4 


6 


2 


6 


3 


-3 


? 


li(ii) 


5 


7 


7 


4 


6 


7 


7 


3 


+1,-2 


8 


8(10) 


5 


8 


7 


8 


6 


8 


8 


3 
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Figure 13: Least recently entered rule 



Step 
Dumber 


Vertex 
(biDary) 


IrectioD 


OotgoiDglrectioiis 
(bold lor cboseii) 


+1 


-1 




-2 


+3 


-3 




-4 


1 


i(iiii) 


i 


1 


i 




i 




i 






2 


1(11) 


2 


1 


















3 


5(11) 


2 


1 







3 




i 




-L+2,|4 


4 


13(1101) 


2 


1 


i 




3 




4 




11-3 


5 


15(1111) 


2 


1 


5 




3 




4 




-14 -4 


6 


11(11) 


2 


1 


5 




3 


6 


1 




-2 




8(11) 


2 


1 


5 




3 


6 


4 




-1 


8 


8(10) 


2 


8 


5 


1 


3 


6 


4 








Figure 14: Least iterations in basis rule 



(binary) 


directioD 


Outgoing directions 
(bold lor cbosen) 


+1 


-1 


+2 


-2 


13 


-3 


+4 


-4 




i 






1 




1 


i 


1 


II +1+3 14 




1 






2 




2 




2 


+2,11+4 


5(il) 


2 




i 


3 


1 


2 




3 


-[m 


13(11) 


3 






4 


2 


2 


1 


3 


+2,-3 


15(1111) 


4 




1 


4 


3 


2 


2 




-1-3-4 


14(11) 


4 




2 


4 


4 


2 


3 




-3 


lO(li) 


4 


3 


3 


4 


4 


3 


4 




+1.-2 


11(1011) 


5 


3 


4 


4 


4 


4 


5 




-2 


9(11) 


5 


3 


4 


5 


4 


5 


6 


3 


-1 


8(10) 


5 


4 


4 


6 


4 


6 


1 


3 





19 



